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ABSTRACT 

This paper presents a study of quasi-steady spherical accretion in the early Universe, before the 
formation of the first stars and galaxies. The main motivation is to derive the basic formulas that 
will be used in a companion paper to calculate the accretion luminosity of primordial black holes and 
their effect on the cosmic ionization history. 

The following cosmological effects are investigated: the coupling of the gas to the CMB photon 
fluid (i.e., Compton drag), Hubble expansion, and the growth of the dark matter halo seeded by the 
gravitational potential of the central point mass. The gas equations of motion are solved assuming 
either a polytropic or an isothermal equation of state. We consider the cases in which the accreting 
object is a point mass or a spherical dark matter halo with power-law density profile, as predicted by 
the theory of "secondary infall" . Analytical solutions for the sonic radius and fitting formulas for the 
accretion rate are provided. 

Different accretion regimes exist depending on the mass of the accreting object. If the black hole 
mass is smaller than 50 — 100 Mq, gas accretion is unaffected by Compton drag. A point mass and 
an extended dark halo of equal mass accrete at the same rate if M ^ 5000 , while smaller mass 
dark halos accrete less efficiently than the equivalent point mass. For masses M x 10^ M©, the 
viscous term due to the Hubble expansion becomes important, and the assumption of quasi-steady 
flow fails. Hence, the steady Bondi solutions transition to the time-dependent self-similar solutions 
for "cold cosmological infall" . 
Subject headings: 



1. INTRODUCTION 

The gravitational collapse of relativistic matter dur- 
ing the early evolution of the universe is thought to 
produce primordial black holes (PBH). Moderately non- 
linear inhomogeneities that become gravitationally un- 
stable after entering the particle horizon have a den- 
sity and radius that is very close to the black hole 
regime. It is easy to show that the Jeans length and 
the Schwarzschild radius of a black hole with mean den- 
sity close to the cosmic value both approach the particle 
horizon radius as the effective sound speed of the gas 
approaches the speed of light. Hence, small density per- 
turbations with 5p/p^Q. 1 — 1 are unstable and collapse 
into PBHs with a mass that roughly equals the mass 
within the par ticle horizon at the redshift of their forma- 
tion (iHawking" 19711: [Carr fc Hawkind ll974l; iMusco et all 
[2005l:THarada fc Carill2no'5f K 

Many physical processes may produce the necessary 
level of inhomogeneity for significant PBH formation. In 
most models PBHs are produced during phase transi- 
tions of the equation of state and, depending on the time 
of their formation, may have masses that range between 
the Planck mass to 10'^ Mq. Although the probability of 
PBH formation is finite, it may be so small to be of little 
or any cosmological interest. PBHs with mass < 10^^ g 
are thought to evaporate emitting Hawking radiation on 
time scales shorter than the age of the Universe. The 
radiation or particles they emit would make them de- 
tectable. Thus observational constraints on the abun- 
dance of evaporating PBHs are rather stringent. The up- 
per limit on the density parameter of PBHs with masses 
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Ig < MpBH < 10^5 g is /3(A£) ~10-20 - lO'^^ at 
the redshift of formation (e.g. ICarill2003l ). On the other 
hand, the existence of non-evaporating PBHs with mass 
> 10^^ g is very poorly constrained. It is conceivable 
that non-evaporating PBHs constitute all or a substan- 
tial fraction of the dark matter in the Universe. The 
existence of PBHs with M ~ 0.1 — 1 M© is of particular 
interest because the MACHO collaboration has detected 
compact objects in this mass range in the Galactic halo, 
const ituting about 20% of the dark matter (jAlcock et all 
|2000| ), althoug h this cla im is weakened by recent results 
(e.g.. [Hamadac he et al.l [2006). PBHs of about 1 M© are 
also of special interest theoretically because they may 
have formed copiously during the quark-hadron phase 
transition due to a temporary so ftening of the cosmic 
equation of state (*Jeda mzi3ll997l ). 

The main motivation for the calculations presented 
this work is to provide the basic formulas for the ac- 
cretion rate onto PBHs that, typically, are thought to 
have masses smaller than 10^ M©. We will show that in 
this mass range the cosmological spherical accretion so- 
lutions are described by the steady equations (i.e., Bondi 
solutions) while for more massive black holes the accre- 
tion equations are time-dependent and the solutions self- 
similar (i.e., secondary infall solutions). 

These formulas will be the foundation fo r further cal- 
culat ions presented in a companion paper (jRicotti et al.l 
l2006l hereafter, Paper II) aimed to understand whether 
the interaction between an, as yet undetected, popula- 
tion of non-evaporating PBHs and the rest of the mat- 
ter in the Universe may produce observable signatures. 
Thus, in this paper we will only focus on the mathemati- 
cal aspect of the spherical accretion problem (i.e., for the 
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mass range we are interested in, the solution of the steady 
"Bondi" problem) without a discussion on the physics of 
the early universe, feedback effects and other complica- 
tions. Nevertheless, in the remaining of this introduction 
we provide some additional background discussion to put 
the mathematical problem into context. 

PBHs can accret e gas and dark matter from a nearly 
uniform medium (|Caril I198H iMiller fc Ostrikeil 120011 : 
iMack et all I2006D well before the collapse of the first 
nonlinear structures (stars and galaxies). Since the den- 
sity inhomogeneities of the gas and dark matter decrease 
with increasing redshift, we expect the angular momen- 
tum of the accreted material (gas and dark matter) to 
be small and hence, the accretion flow to become in- 
creasingly spherical. At redshifts greater than z ^ 100, 
CMB photons are strongly coupled to the gas because of 
Compton scattering with the residual electrons left over 
after the epoch of recombination [xe — 10~*). The gas 
behaves as a viscous fluid and the accretion onto PBHs 
can be substantially reduced. The photon viscosity may 
also be impo rtant in rem oving angular momentum from 
the gas ('e.q.. lLoebiri993D . We will show that the Hubble 
expansion also acts as a viscous term that tends to reduce 
the accretion rate, independent of the ionization fraction 
of the gas. Finally, if massive PBHs do not constitute the 
bulk of dark matter, their gravitational potential seeds 
the formation of a dark halo that is expected to grow 
over time becomin g ^ 100 times mo re massive than the 
PBH in its center (jMack et al.ll2006[ ). 

The gas accretion onto PBHs produces ionizing radi- 
ation and an X-ray background that may be sufficiently 
large to heat the cosmic gas and keep it partially ionized 
after recombination. Thus, PBHs may imprint signa- 
tures on CMBR w hich are incompatib le with observa- 
tions by WMAP3 ijSpergel et al.ll2006D and such obser- 
vations can be used to constrain PBH masses and abun- 
dances. Several ingredients need to be considered when 
modeling ionization by accreting PBHs, including their 
proper motion and feedback effects. The X-rays emitted 
by accreting PBHs, for example, increase the tempera- 
ture and fractional ionization of the gas and are able to 
reduce the accretion rate. Thus the global accretion rate 
onto PBHs and the ionization state of the cosmic gas 
are coupled and need to be calculated self-consistently. 
Such calculations will be presented in a companion paper 
(Paper II). 

The present study focuses on the derivation of the ba- 
sic analytical relationships for the gas accretion rate as a 
function of redshift, temperature, density and ionization 
fraction of the gas. We find solutions for the spherical 
accretion problem including the following cosmological 
effects: i) the effect of the viscous terms due to Comp- 
ton drag and Hubble expansion; ii) we consider both the 
case of an isolated PBH and the case in which the PBH is 
"clothed" by a dark matter halo with power-law density 
profile. We discuss the relationships between the steady 
Bondi type solutions and the time-dependent-self-similar 
solutions of secondary infall. Analytical expressions for 
the accretion rate as a function of the "viscosity parame- 
ter" are provided for either a "naked" or "clothed" PBH. 
We also consider the different cases of isothermal equa- 
tion of state or polytropic equation of state with index 7 
f or the gas. 

lUmemura fc Fuku3 ()1994l ) have previously investi- 



gated the effect of radiation drag on spherical accretion 
considering three cases for the gravitational potential: a 
central point mass, a dark halo with constant potential 
and a self-gravitating gas sphere. Their work is similar 
to the present work for the case of accretion onto a point 
mass. Our work differs from the previous ones in that it 
addresses the growth of a dark matter halo around PBHs, 
it discusses the validity of the Bondi accretion regime 
and the transition to time-dependent solutions. Finally, 
and most importantly, we derive analytic expressions for 
the accretion eigenvalues, which are a key ingredient to 
model the cosmic ionizatio n by PBHs presented in Pa- 
per II. lTsuribe et al.l (|1995l ) has also investigated a spher- 
ically symmetric accretion flow with an external radia- 
tion drag onto black holes but for masses that are larger 
(> 10*"^ M0) than the ones discussed in this paper. We 
will show in § 12.11 that in this mass range the steady- 
state assumption that defines the "Bondi problem" fails 
and the solutions transition to a class of time-dependent 
self-similar solutions. 

Another interesting mechanism that may lead to the 
formation in the e arly un i verse o f massive black h o les ha s 
been proposed bv iLoebl (|1993f l: lUmemura et al.l (|1993f l. 
These works discuss in detail the effect of CMBR in the 
accretion process during a quasi-spherical collapse of a 
self-gravitating gas clouds and the consequent formation 
of massive black hole in the early universe, after recom- 
bination. In this scenario the removal of angular mo- 
mentum due to Compton drag in collapsing rare density 
perturbations at redshift z ^ 400 — 1000 leads to the 
formation of massive black holes with masses M ^ 10^ 
Mq. Radiation feedback effects associated with the ac- 
cretion process are also considered. These works differ 
from the present one in several ways. They focus on the 
formation of black holes rather than accretion onto a pre- 
existing black holes and the typical masses of the black 
holes formed by this process are much more massive than 
the ones considered in the present study. The accretion 
is not steady and no Bondi solution exists. Although the 
black holes formed in this scenario should be very rare, 
;Sasaki & Umemura (199f|) have estimated that accretion 
onto such black holes would produce enough radiation to 
reionize the universe \iy z '-^ 150. 

Several authors have investigated the effect that 
Compton drag and Compton heatin g has on g as accre- 
tion in many different co nt exts ( e.g., Fukue fc U memur a 
T99I iMi ncshiec et all [19981 iCiotti et all l200l 
Park fc Ostriker 200l . Recently, IWang et all (|2006l ) 
presented a discussion of feedback in different accretion 
regimes, including a Bondi sphere in the case of a 
black hole accretion at high redshift. They point out 
that Compton heating can rapidly halt supercritical 
accretion onto remnant black holes of 10'^ Modot from 
Population HI stars. Such feedback would prevent seed 
black holes to grow to masses typical of supermassive 
black holes on timescale shorter than a Gyr. Thus, 
they conclude that accretion onto seed black holes from 
population III stars is not a viable mechanism to explain 
the existence of quasars at redshift z ^ ^. In the present 
paper the accretion onto PBHs is highly subcritical 
and Compton heating is unlikely important. However, 
feedback processes and other physical processes will be 
addressed in detail in Paper II. 

The organization of this paper is as follows. In § [2] 
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we introduce the basic equations for accretion around 
a point mass in comoving coordinates. In § |3] we solve 
numerically the accretion equations and provide analyt- 
ical fits for the accretion eigenvalues. In § |4] we repeat 
the calculations assuming the gravitational potential of 
a dark halo with power-law density profile. A summary 
and discussion is presented in § [5] 

2. BASIC EQUATIONS 

At high redshift, repeated Compton scattering between 
free electrons and CMB photons - assumed here to be ho- 
mogeneous and isotropic - act on the proper velocity of 
the gas, V, as a viscous term. The forces exerted by the 
CMB photons are opposed to the motion of the fluid ele- 
ment with respect to the phot ons' rest frame , producing 
a mean restoring acceleration ()Peebleslll98Cl[ ) 



comoving radius, x, 



dv 



4 XeCrrU, 



3 



cmb 



V — —I3v, 



where Xe is the electron fraction, ctt is the Thompson 
cross section to electron scattering and Ucmb = ^T^mb 
is the CMB energy density. By expressing the CMB 
temperature as a function of the redshift, z, we have 
(3 = 2.06 X 10-23a;e(l + z)^ s"^. Note that the typical 
time scale for Compton drag is nip/me ^ 1000 times 
longer than the timescale for Compton cooling/heating. 

The "Bondi" problem is described by the stationary 
equations of mass and momentum conservation (includ- 
ing the pressure, gravitational potential and viscosity 
term) and the equation of state of the gas. In spheri- 
cal coordinates the equations are: 



M, 



Anr pv, 



dv _ 1 dP 
dr p dr 

P = i^pT. 



GM(<r) 



I3{z)v, 



(1) 



Here, v, r and p are the velocity, distance from the cen- 
ter of the gravitational potential and gas density. The 
"dot" represents the time derivative, thus M is the mass 
accretion rate. The gas pressure, P is assumed to be a 
power law of the density with exponent 1 < 7 < 5/3. 
This equation, know as the polytropic equation of state, 
substitutes the energy conservation equation. The case 
7 = 1 corresponds to an isothermal system in which the 
cooling time is much shorter than the dynamical time 
and the opposite limit, an adiabatic system, is obtained 
for 7 = 5/3. 

So far we have neglected the cosmological terms due 
to the Hubble expansion and dark matter. Also, we have 
implicitly assumed that the redshift-dependent viscous 
term, (3{z), changes slowly with respect to the time scale 
necessary to achieve stationarity: tcr <^ tn, where tcr 
is approximately the sound crossing time and tn is the 
Hubble time. We will return to this point in the next 
section. 

2.1. Hubble expansion: transition to self-similar 
solutions 

In order to include the cosmological term due to the 
Hubble expansion we use the standard procedure of writ- 
ing the hydrodynamic equations ([T|) in a comoving frame 
of reference. This is done by expressing the radial co- 
ordinate, r, and the radial velocity, v, in terms of the 



and peculiar velocity, Vpi 

a{t)x, 

dt^ 



^-,{ax) 



Hr 



(2) 



Here, a{t) is the scale factor, H = a/a is the Hubble 
parameter and the peculiar velocity is defined as Vp — ax. 
The stationary {06 /dt — 0, dvp/dt = 0) equations for the 
overdensity, 6 = p/j), and peculiar velocity, Vp, are 



Mr, 



Attx'^{1 + S)vp 



v^Ov^ ^_ldP_ G\AMi<a.)] _^^ j^. 
a ax p adx [ax)'^ yt-' ' j p 



(3) 



Multiplying both sides of the mass conservation equation 
(top equation in [5]) by 'p{t)a{t)'^, we recover the same 
equation as in ([1]), where Mg — 'p{t)a{t)^Mp is now a 
function of time (or redshift). We also notice that 8/ dr = 
a-^d/dx and that AM(< r) = M(< r) - ^ifpr^ jZ « 
M(< r). Thus, neglecting the self-gravity of the accreted 
gas, the momentum conservation equation (the bottom 
equation in [3] ) has the same form as in ([T]) after replacing 
/? with the effective viscosity jS'^^^ = {[3 + H): 



Mg 

V ^ 

^ dr 

P 



= 4'Kr^p{z)vs, 

GM(<r) _^e//(^) 

= Kp'^. 



(4) 



If we neglect the gas self-gravity, the cosmological ex- 
pansion does not change the basic form of the equations 
but introduces: (i) an additional viscosity term in the 
momentum equation and (ii) the boundary conditions 
at r ^ cxD for the density Poo{z) and the sound speed 
Cs.oo{z) are a function of redshift. Although we are as- 
suming stationarity in the comoving frame of reference, 
the flow is time-dependent in physical coordinates. The 
assumption of steady flow is justifled as long as the time 
scale to achieve the stationary solution in the comoving 
frame of reference is shorter than the Hubble time. We 
have already made this assumption when we have con- 
sidered the Compton viscosity, I3{z). 

2.1.1. Stationarity 

The typical scale of the system is the Bondi radius 7'b, 
hence, the sound crossing time is tcr ~ h/cs^oo- If tcr < 
tn the Bondi solution is valid because the cosmological 
terms vary slowly with respect to the typical time scale 
of the system and a quasi-steady solution can be found. 

The analysis of the dimensionless accretion equations 
(see § [3]) show that the viscosity term due to the Hubble 
expansion becomes important when Hrb/cs^oo > 1- This 
condition is equivalent to the constrain tcr > tn- Thus, 
the assumption of stationarity, implicit in the derivation 
of the Bondi solution, fails when the Hubble term be- 
comes important. In §[5l assuming a static PBH and ne- 
glecting feedback effects, we estimate that the hypothesis 
of stationarity fails for PBH masses Mpbh ~ 2 x 10"* M©. 
This is a conservative estimate because if we consider the 
proper motion of PBHs and feedback effects we derive 
a larger critical mass. However, the existence of PBHs 
more massive than 10'' — 10^ M© is unlikely, so the cases 
for which the Bondi solution does not apply are not par- 
ticularly relevant for the aim of the present study. 

2.1.2. Gas Self-gravity 
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Fig. 1. — (Left). Dimensionless accretion rate A onto a point mass as a function of the effective viscosity, /3, produced by Compton drag 
and Hubble expansion. We assume an isothermal equation of state for the gas. The points show the result of numerical integration and the 
curve shows the fit with the relationship in equation (|9j. {Right). The same as in the left panel but assuming a polytropic equation of state 
for the gas with 1 < 7 < 5/3. Instead of the accretion eigenvalue A, here we show the sound speed at the sonic radius, Cs.cr (normalized 
to its value neglecting viscosity, = 0), as a function of the effective viscosity 0. The accretion rate can be calculated from equation (|8j. 
The curves from top to bottom refer to 7 = 1.1, 1.2, 1.3, 1.4, 1.5, 1.6 and 1.65, respectively. 



The assumption of neglecting the self-gravity of the 
gas is justified if M ^ Mgas, where Mgas is the total 
gas mass contained in the halo around the black hole. 
It is easy to estimate the gas mass, M^^^, in the limit 
in which the viscosity terms are negligible {i.e., for the 
classic "Bondi" solution). This estimate is a lower limit 
for M because M > M^^^ > Mgas- By integrating 
the density profile for the Bondi solution we find that 
^^gas ^ 87r/3poo'"b, whcrc rfc — GM/c^^^ is the Bondi 
radius. Thus, for M <C M'^'^ = (SttG^ /S)-^/'^cl^pZo^^^ 
we can neglect the gas self-gravity. In the case of PBHs 
accreting from the IGM between the redshift of recombi- 



nation and decoupling (100 , 



1000) we have Cs 



5.74km s-i[(l -t- z)/1000]i/2 and p^o = PcAil + z) 



Assuming Q,i,h^ — 
later times (30 < 



0.02 we find M^'' - 3 x 10^ M©. At 
z ^ 100) the critical PBH mass may 
become larger than M'^^ because the temperature of the 
IGM may increase due to the radiation emitted by the 
gas accreting on the PBHs. The condition for steady ac- 
cretion is more restrictive than the one found here. Thus, 
if steady solutions exist, the gas self-gravity can always 
be neglected. 

The solution of the Bondi problem including the gas 
self-gravity has been studied in the context of stellar 
structure. The solution is known to be unstable for 
7 < 4/3. In this limit we would find that all the gas 
in the universe is eventually accreted onto PBHs. In the 
regime discussed here only a small fraction of the cosmic 
gas is accreted. 



the growth of stru ctures seeded by the presence of PBHs 
(jMack et al.ll2006f) . The total energy of a representative 
patch of the Universe at critical density is zero. Thus, 
the volume is marginally stable and the introduction of 
any perturbation, such as a PBH, produces the collapse 
o f the region. 

iMack et al.l (|2006[ ) found that a dark matter halo about 
100 times more massive than the primordial black hole 
at its center is accreted by redshift z ^ 30. The accretion 
takes place mostly after the epoch of "matter-radiation" 
equivalence. Calculations of the gas accretion rate into 
the gravitational potential produced by the dark matter 
halo will be considered in § |4l When a dark halo grows, 
the gas self-gravity can be neglected with respect to the 
dark matter potential because the total gas mass is about 
five times smaller {i.e., the cosmic value) than the dark 
matter. 



3. SPHERICAL ACCRETION ONTO A POINT MASS 
3.1. Sonic point 
It is convenient to write equation ^ in dimensionless 



r/n, M = Vs/Cs, Cs = Cs/Cs, 00, P = P/Poo, 

s^oo and Poo are the sound speed and gas density 



units, X 
where c 

at radii r —^ 00, respectively. The scaling constant for 
the radius, viscosity and accretion rate are: 



GM, 



n 



bh 



n 



A = 



M„ 



(5) 
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2.1.3. Dark matter self-gravity 

The assumption of neglecting self-gravity of the ac- 
creting material is a good approximation for the baryons 
but not for the coUisionless dark matter component. Un- 
less all dark matter is constituted of primordial black 
holes {i.e., Qdm — ^bh) we need to take into account 



We obtain the familiar set of dimensionless equations: 



A = p^~Mx^ 



M 
M 



M'-'-l 



I3M 



(6) 
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The transonic solution crosses the sonic point {M. = — 1) 
at the critical radius: 



+ \ 1 + I3c. 



3 
s.cr 



I3c. 



3 

■s,cr 



(7) 



where Cs^cr is the dimensionless sound speed at the crit- 
ical radius. In the limit of small effective viscosity (i.e., 
/3^c^,cr) find Xcr (2Cs cr)~^- the limit of large 
effective viscosity we find Xcr (cs,c/3)~^^^- The eigen- 
value, A, for the transonic solution is 



Pc 



7 + 1 



(8) 



where pcr is the dimensionless density at the critical ra- 
dius. 

3.2. Accretion eigenvalues 
3.2.1. Isothermal equation of state 

Since we did not find analytic solutions for equa- 
tion ([6]), we integrated the equations numerically. We 
found the asymptotic solutions of the equations in the 
limits of negligible viscosity, /3 — > 0, and large viscosity, 

— > oo. In the limit of small viscosity we find the classic 
Bondi solution A — > exp(3/2)x^^, where Xcr — 0.5. In the 
limit of large effective viscosity, the asymptotic analysis 
of equation ([6]) shows that pcr 1 and that A — > /3~^. 
We have used the asymptotic solutions for A to guess the 
form of the appropriate fitting function of the numerical 
results. The function 



A — exp 



9/2 



3 + [3 0-75 



(9) 



gives a good fit for the numerical results as shown in 
Fig.mieft). 

3.2.2. Polytropic equation of state 

The asymptotic solutions for the cases of negligible 
effective viscosity and large effective viscosity can eas- 
ily be found. Neglecting the viscosity terms we have 
Cs,cr — > [2/(5 — 37)]^/^. In the limit of large viscosity the 
asymptotic analysis of equation ^ shows that pcr 1, 
Csxr 1 and A — > /3^^. We solve equation ([6]) numeri- 
cally and find a parametric fit for the numerical results, 
with the correct asymptotic behavior for the sound speed 
at the critical radius. We also make sure that in the limit 
7 — > 1, the parametric fit for Cg cr and A have the same 
expression as that found for the isothermal case. We find 



Cs,cr(/^) " 



3 + /3« 



^ = 0.75 



7 + 1 



5-37 



(10) 



where Cs^criO) ~ [2/(5 — 37)]^/^ is the sound speed at the 
critical radius for /3 = 0. The fitting curves (solid lines) 
and the numerical integrations (shown by the symbols) 
are shown in Fig. [TJright). The eigenvalue A can be cal- 
culated for any given value of the polytropic index, 7, 
combining equations ([8]) and (flO)) . 



4. SPHERICAL ACCRETION INTO A DARK MATTER 
HALO 

If dark matter is entirely composed of PBHs each of 
them would only accrete gas. In all the other cases a dark 
halo of WIMPs forms around each PBH. In this section 
we examine the gas accretion rate into the gravitational 
potential of a spherical dark matter halo with power- 
law density profile. The assumption of a single power- 
law for the density profile is partially dictated by the 
need for simplicity but it is also justified theoretically as 
elucidated below. 

PBHs accumulate a spherical halo of dark matter 
that grows in mass proportionally to (1 -|- z)~^. Non- 
negligible growth of the dark matter halo begins at 
matter-radiation equality, and by redshift z ~ 30 the 
black hole is embedded in a h alo about 100 times its 
original mass (|Mack et al.ll200l . The numerical values 
of the halo mass and the comoving turnaround radius 
are: 

M. = 3Aw(i±i) \ (11) 



rta., 



:58 pc 



1 Mr. 



1/3 



(12) 



The theory of secondary infall predicts that the dark 
halo has a radial density profile that is well approxi- 
mated by a single power-law cusp and a truncation of 
the density profile at a fraction of the turnaround radius. 
All of the halo mass is contained within the cusp. Two 
types of time-dependent, self-similar solutions have been 
found for the dark rnatter profile ([G unn & Gottl 119721 : 
iFillmore fc Goldreichlll984l : lBertschiVer.l985i ): the tran- 
sition from one type to the other depends on the amount 
of angular momentum of the infalling matter: 

a) If the flow is highly radial such that the PBH at the 
center absorbs all of the infalling matter, then almost 
all the mass in the accreting halo is eventually incorpo- 
rated into the central black hole. This infall solution is 
characterized by the free-fall boundary condition at the 
center and pressureless matter at infinity. The solution 
is also valid for "cold" gas accretion because there is no 
shell crossing and no standing shock forms. The self- 
similarity of the solution is valid if Mh > Mpbh, when 
the matter is self-gravitating. Near the center the halo 
mass is negligible compared to the black hole mass, so the 
matter is free falling and the density profile is p cx r~^^'^. 
This slope is the same as for the Bondi solution inside 
the sonic radius. In the outer parts the solution is very 
different from the Bondi solution. The radial density pro- 
file, due to the Hubble expansion, decreases with radius 
propotionally to r~^. The radius rh of the "free-fall" 
cusp depends on the details of the collapse but is typi- 
cally rh/rta < 1 (|Mack et al.ll2006D . 

b) If the fiow is quasi spherical but the angular momen- 
tum of the accreting dark matter prevents direct accre- 
tion onto the central black hole, the particles pass near 
the center without being absorbed. For the case in which 
the accreting matter is collisional, a standing shock inside 
the turnaround radius stops the gas infall within the core. 
Remarkably, the solution for coUisionlcss matter is very 
similar to the 7 = 5/3 solution for gas accretion. The 
density profile has slope a — 2.25 a,t r < rh = 0.339rtQ. 
For r > rh the overdensity is close to one. 
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Fig. 2. — (Left). Dimensionless accretion rate A into a spherical dark matter halo as a function of the cfEcctive viscosity, 0, produced 
by Compton drag and Hubble expansion. Each curve, from bottom to top, refers to halos of increasing radius. If the radius of the halo is 
smaller than the accretion radius, r^, defined in the text, the accretion rate into an extended halo equals the one calculated for a point 
mass. We assume isothermal equation of state for the gas and a power-law density profile of the dark halo with slope a = 2.8. From bot tom 
to top the symbols show the result of numerical integration and the curves show the fits with the parametric relationship in equation (I18I I 



for values of the dimensionless halo radius rg/vh = X 



: 2, 1, 0.5, 0.25 and 0.016. (Right). The same but for a = 2.4. 



Estimates of the angular momentum of the accreting 
dark ma tter onto PBHs a re presented in a companion 
paper (iR icotti et al.l 120061 ). Based on the results of the 
angular momentum calculations case (b) is the most rel- 
evant. Nevertheless, if a fraction of the dark matter is 
captured by the central PBH, we expect that the den- 
sity profile has properties intermediate between case (a) 
and (b) . The truncation radius would be smaller than 
Tta/i, and the effective slope of the inner density profile 
steeper than 2.25, in order to take into account the non- 
negligible mass of the central black hole. Indeed a halo 
with slope of the density profile a = 3 has a potential 
that approaches the one of a point mass with M ~ Mh- 
For the sake of generality we take this effects into ac- 
count by treating the core radius Vh and slope of the 
core density profile as free parameters. 

We calculate of the gas accretion rate by replac- 
ing the gravitational acceleration for a point mass 
{gbh = —GMbfi/r'^) in equation ([6]) with that for an ex- 
tended dark halo with power-law density profile: p = 
^Pdm{z)[r /rh]~°' ■ Here, pdm{z) is the mean dark matter 
density at redshift z and A is the overdensity of the halo 
&i r = Th- We assume that all the mass of the halo is 
within the radius rh. Thus, the gravitational acceleration 
is 

gdmir) = —ram\{r/rhY, 1], (13) 

where p = 3 — a. Note that the halo radius is smaller 
than the turnaround radius: r/j < rta- We explore the 
range of values 2 < a < 3 for the log-slope of the density 
profile (i.e., < p < 1). If p ^ the gravitational ac- 
celeration approaches the point-mass value. If the halo 
radius is Th < Xcr = 0.5r(, = 0.5GAf/i/Cg the accretion 
eigenvalue is given by equation ^ , the same as replacing 
the halo with a point-mass M = Mh- Thus, it is appar- 
ent that the accretion rate onto an extended halo is a 
function of the dimensionless parameter x — Tb/^h and 



the parameter p. The case p = 1 describes the density 
profile for an isothermal sphere for which the circular ve- 
locity is a constant. This case is singular because the 
Bondi radius, for which Vdri^b) = Cs, cannot be defined. 

4.1. Sonic point 
After substituting g{r) in equation ([4]) with the expres- 



sion given in equation (jl3p for gdmif), we rewrite the 
equations in dimensionless units, x = r/rs, M = Vg/cs, 
Cs = Cs/cs,oo, P = p/poa, where c^^oo and poo are the 
sound speed and gas density respectively, at radii r — > oo. 
Note that the scale radius, tb, for spherical accretion 
onto an extended halo of mass Mh is not defined analo- 
gously to point-mass Bondi radius, = GM^j c?^^^, of 
equivalent total mass. It is convenient to express the di- 
mensional constants for an extended halo as a function 
of the parameter x — Tbl^h'- 



p)^^ < oo and < p/(l 



(14) 



where 1 < (1 — p)^^ < oo and < p/(l — p) < oo for 
2 < a < 3. The dimensionless viscosity and accretion 
rate are defined as: 



-X 



In these units we obtain the dimensionless equation 



(15) 



M 
M 



(7-1) 



2± 
2 



1 ( 1 



(16) 



- 1 

The transonic solution crosses the sonic point {A4 = — 1) 
at the critical radius Xcr,dm, found by solving the equa- 
tion 

where Cs^cr is the dimensionless sound speed at the crit- 
ical radius. It is useful to find the asymptotic behav- 
iors of the critical radius in the limit of zero viscosity. 
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Xcr,dm ^ i'^^t cr) ^^^^ ^^d in the limit of large vis- 
cosity Xcr,dm {cs,cr$)^^^''^^^'' ■ Tliesc asymptotic be- 
haviors will guide us to find the functional form of the 
fitting formula. 

4.2. Accretion eigenvalues 

For the sake of simplicity, given that we have intro- 
duced two new free parameters that describe the halo 
potential, we will only discuss the case of an isothermal 
equation of state for the gas. Physically, the isothermal 
equation of state is the most relevant case because the 
CMB photons that exert Compton drag are about 1000 
times more efficient in keeping the gas temperature con- 
stant, and close to the CMB value. 

4.2.1. Isothermal equation of state 

The accretion eigenvalue can be fitted with the follow- 
ing approximate formula 



A" 



9/2 



l5 + (3 



0.75 



^ cr.dm 1 



where the functions 
fx- 



1 + 1.25/3^, 
:[exp(2-x)]" 



(18) 

(19) 
(20) 



are small corrections to A that vanish for small values of 
p and X — > 2. This formula cannot be used without an 
expression for Xcr,dm- The function fp has been derived 
from the asymptotic analysis of the equation for the crit- 
ical radius Xcr,dm- Assuming an isothermal equation of 
state, the sonic radius approaches Xcr,dm 2^^/'-^^^^ in 
the limit of negligible effective viscosity, and Xcr,d7n 
/j-i/(2-p) ^ijg limit of large effective viscosity. Thus, 
the sonic radius of the transonic solution for accretion 
into a halo potential is related to that of a point-mass 
(with M = Mfi) by the approximate fitting formula 



Xcr,dm 



Xc 



1/2 



where 



5/3 



1 + 10(3 



lO(l-p) 



(21) 



(22) 



Combining equations (PT|) and we can relate the 
accretion eigenvalue for an extended dark matter halo to 
the eigenvalue in equation ([9]) for a point mass AI — M^'- 



A^ 
A 



(23) 



The fits using the parametric function equation (|18p to 
the numerical integrations are shown in Fig. [5] for differ- 
ent values of a and the parameter x- 

It is important to note that the dimensional units for 
the accretion rate, radius and effective viscosity are dif- 
ferent for the cases of a dark matter potential and a 
point-mass. Thus the dimensional accretion rate into a 
halo is reduced by a factor x^^^ (i~p) and the effective vis- 
cosity increases by a factor x^^^^^^-* with respect to the 
point-mass case. To summarize, the accretion rate into 



a halo of mass Mh is related to the accretion rate onto a 
black hole of equivalent mass Mbh by the relationship 



A'': 



„h — 



TT^A^'\ 



. 1/2 , 



where 



T = 1 + 10/3^ 



exp (2 - x) ( 2 



(24) 
(25) 

(26) 



(27) 



Note that and not should be used in the analytic 
expressions for A'''' and rcr- 

5. SUMMARY AND DISCUSSION 

This paper examines the classic Bondi problem of 
steady, spherical gas accretion in a cosmological frame- 
work. The main motivation of the present work is to de- 
termine the accretion rate onto primordial black holes in 
the early universe, before the formation of the first galax- 
ies and stars. We estimate the importance of cosmolog- 
ical effects such as the gas viscosity due to the coupling 
of the gas to the CMB photon fluid (Compton drag), 
the Hubble expansion and the formation of a dark halo 
around the black hole. In a companion paper we use the 
results of the present work to model cosmic ionization by 
the X-rays emitted from accreting PBHs. Comparisons 
between theoretical models and WMAP3 d ata allow us to 
const rain mass and abundances of PBHs (jRicotti et al.l 
[2006h . 

The cosmological Bondi problem is tightly re- 
lat ed to the the o ry of "secondary infall" develope d 
bv |Gunn fc Got^ g97l) : iFillmore fc GoldreichI ([1981 : 
iBertsc hingcr ( 1985() . In the standard cosmological model 
the baryons are a small fraction of the dark matter, hence 
the evolution of dark matter is unaffected by the presence 
of the gas. We use the results of the secondary infall the- 
ory to determine the shape of the gravitational potential 
into which the gas is infalling. In addition, when the vis- 
cosity term due to cosmic expansion becomes important, 
the assumption of steady flow for the gas, implicit in de- 
riving the Bondi solution, becomes invalid. In this limit, 
time-dependent infall solutions should be used instead. 
In § 15.31 we discuss the relationships between the Bondi 
and the secondary infall solutions and the conditions of 
validity for each of them. 

5.1. Accretion onto a point mass 

The physical values for the accretion rate and the crit- 
ical radius are easily calculated from the dimensionless 
ones. The values of the Bondi-Hoyle radius, rf,, and ac- 
cretion rate are 

/ M 



«43.3 pc 



1km s-i 



/ M 



1km s" 



(28) 
t)(29) 



M«A(3.7xlO-gs-),.,..^^^,^^ 

respectively (|Bondi fc Hovld[l944[ ). The accretion rate 
eigenvalue. A, is a function of the dimensionless "effective 
viscosity", (3 — /3'^^-^tcr, where the Bondi radius crossing 
time is tcr — 7'h/cs,oo and (3'^-^^ = I3{z) + H{z). In the 
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Fig. 3. — (Left). Dimensionless accretion rate A into a spherical dark matter halo (solid curve) and a point mass (dashed curve) as a 
function of the mass of the object (PBH or dark halo) at 2 + 1 = 1000. We have considered only the viscosity from the Hubble expansion 
and a halo with power-law slope a = 2.8. (Right). The same as the left panel but considering only the viscosity due to Compton drag. 
The different curves refer to gas with fractional ionization Xe = 1, 0.1, 0.01 and 0.001 as shown by the labels. 



relevant redshift range {i.e., before z ~ 30) the Compton 
drag term is (3{z) oc Xe (1 + z)4 and the Hubble param- 
eter is H{z) = 3.24 X 10-i8s-i(17 /t2)i/2(l-Hz] 3/2^ 

follows, assuming il^/i^ = 0.127 (ISpergel et al.ll2006D . 

that P'^ff = H{z){l + 1.78a;e[(l + z)/100]^/^} and the 
Compton term is dominant over the Hubble term at red- 
shifts 



Assuming that the gas temperature equals the CMB 
temperature, the viscosity term due to Compton drag is 
negligible for 

1000 ) 



M, 



70 Mf^ 



pbh 



The gas viscosity can be neglected if (3^^ Her < 1, where 



(0.224 Myr) 



M 
10* M, 



.5.74 km s" 



Including the Hubble flow and Compton drag terms, the 
expression for the dimensionless effective viscosity is 



^ ( 10* Mq ) (looo) (5. 



74kms" 



0-257+ 1.45 



5/2' 



(31) 



Let us consider first the viscosity term due to the Hub- 
ble expansion. We can neglect the Hubble flow when 
tcr < tn — H{z)~^ or, equivalently, when vh/cs^oc < 1, 
where vh = f'}jH{z) is the typical Hubble velocity of 
the gas at the Bondi radius. If we assume that the 
gas temperature equals the CMB temperature, c^^oo = 
1.1 km/s(r„„6/100X)i/2, we find that t^r < tn for 



Mpbh < 2.7 X lO'' M0. 



(32) 



Hence, the accretion onto primordial black holes less 
massive than 3 x 10^ M© (or more massive if the gas 
temperature is Tgas > Tcmb due to photo-heating from 
accreting PBHs) is weakly affected by the Hubble ex- 
pansion. But, as discussed in detail in § 15.31 when the 
Hubble viscosity becomes important the hypothesis of 
quasi-steady flow fails and the Bondi solution transitions 
to the time-dependent secondary infall solution. 



Before recombination, PBHs with mass Mpbh ^ 50 — 100 
Mq are unaffected by Compton drag. Instead, the ac- 
cretion rate onto PBHs more massive than 100 — 1000 
Mq is negligible before recombination. In the standard 
cosmological model the residual ionization fraction af- 
ter recombination is ^ 10"'*, hence, after recombina- 
tion, Compton drag is negligible for any reasonable PBH 
mass. Only if a sufficiently large number of PBHs or 
other sources of ionizing radiation increase the fractional 
ionization above x^ ^ 1.47 x 10~^ the Compton viscosity 
becomes important (i.e., dominant with respect to the 
viscosity term due to the Hubble fiow). 

The eigenvalue A can be derived from the equation of 
mass conservation ([8]) that depends on the polytropic 
index, 7, the sonic point, Xcr [given by equation ([7])], and 
the sound speed at the sonic point, Cs,cr- We provide 
fits to the sound speed at the sonic point (see Fig. [1]) in 
equation (jlOp . In the limit 7 = 1 (isothermal equation of 
state) the accretion eigenvalue is given by equation ^ 
and depends only on Xcr. 

5.2. Accretion into a dark matter halo 

If the radius of the dark matter halo, r^, is within the 
Bondi critical radius, the accretion rate is unchanged 
with respect to the case of a point mass. Quantita- 
tively, the accretion rate is reduced approximately by 
(x/2)2p/(i-P) with respect to the point mass case where 



2^ = ZL = 3.3x10-3 
2 2rh 



1 



1000 



5.74kms-iy 



.5.74kms-i- 
(33) 

If x/2 > 1 the extended dark halo accretes as a point- 



mass with AI — Mh. Assuming Tgas = 



crab 1 



which 
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is a good approximation for z > Zdec ~ 100, we have 
approximately 

X ( Mh\i 

where = 5196 Mq. 

In Fig. we compare the dimensionless accretion rate 
at z = 1000 for a point mass (dashed curves) and an 
extended halo (solid curves) as a function of their mass. 
The left panel shows the case in which the effective gas 
viscosity is dominated by the Hubble flow and the right 
panel shows the opposite case in which Compton drag 
dominates. The results for an extended dark halo shows 
that the accretion rate peaks at M ^h"- This can 
be easily understood due to two competing effects. The 
gas effective viscosity becomes larger with the increasing 
mass of the accreting object, thus the accretion eigen- 
value A is reduced as the mass increases. However, the 
accretion rate eigenvalue for an extended dark matter 
halo is increasingly reduced with respect to the point 
mass case as the halo mass decreases. If Mf^ lO'' Mq 
we have that > r^Q. We will show in i) l5.3l that in this 
regime the Bondi solution becomes invalid and the gas 
flow transitions to the time-dependent self-similar solu- 
tion described by the theory of secondary infall. 

The ratio between the accretion rate onto a PBH 
clothed with a dark halo and the same PBH without 
dark halo ( "naked" ) is approximately 

4^ = Max I (^) ^ [1 (l±l)] , l| , 

and is very sensitive to the slope of the halo den- 
sity profile. At a given redshift z, the effect of the 
halo growth becomes important in increasing the ac- 
cretion rate if PBHs have a mass Mpbh > M^^Kz + 
1/1OOO)/0J^'^~''-'/'^^. For example, assuming a steep 
dark matter density profile with p = 0.2, any PBH 
"clothed" in its dark halo, with a mass Mpbh > 
2.37 M0[(3/(;ii)(l-|-z)/lOOO]'^ accretes significantly more 
and its accretion rate increases rapidly with time (cx 
(1 + z)^^/^) with respect to a "naked" PBH of the 
same mass. As noted in § HI density profiles steeper 
than a = 2.25 may be used to describe a halo with 
Mfi ^ Mpbh, in order to take into account the poten- 
tial of the central PBH. Thus, at redshifts z ^ 500, when 
Mfi ^ Mpbh , the fiducial value of the density profile slope 
is a = 2.25 according to the secondary infall theory. In 
this case {p = 0.75 and (pi = 3) the effect of the dark mat- 
ter halo becomes significant for relative massive PBHs: 
Mpbh > 1000 Mo[(l + z)/1000]l■^ or Mpbh ^ 32 Mq at 
redshift z ~ 100. 

5.3. Transition to time- dependent secondary infall 
solutions 

The cosmological Bondi solution rests on the assump- 
tion of steady and non-self-gravitating spherical gas ac- 
cretion. In comoving coordinates the gas is at rest at 
infinity with sound speed, c^^oo- Conversely, secondary 
infall solutions describe self-similar and time-dependent 
spherical accretion in an expanding universe. They de- 
pend on the assumption that the matter is self gravitat- 
ing and its sound speed at infinity is c^^oo ("cold" 



cosmological accretion). Thus, the regime of validity of 
secondary infall solutions is complementary to the Bondi 
solutions. The accretion is effectively "cold" if the sound 
speed at infinity is much smaller than the Hubble flow: 
Cs,oo ^ TbH. This condition is equivalent to tcr < tn] 
that is, when the hypothesis of quasi-steady flow fails for 
the Bondi problem. As shown in § 15.11 for an accreting 
point mass the transition from steady to self-similar so- 
lutions is for Mpbh ^ 2 X 10^ Mq. This case is of little 
interest because the formation of PBHs with such a large 
mass is unlikely. Moreover, we show later that the accre- 
tion rate is super-Eddington, thus regulated by complex 
feedbacks, the modeling of which is beyond the scope of 
this paper. 

When PBHs do not constitute the bulk of the dark 
matter a dark halo accumulates around them. Assum- 
ing, for simplicity, Tgas — Tcmb we find that if Mh <L 10^ 
Mq the Bondi radius is larger than the dark matter 
turnaround radius: rp > rta- This result is clearly incor- 
rect because the gas cannot be accreted from radii larger 
than the dark matter turnaround radius and is simply 
a consequence of neglecting the Hubble expansion. This 
suggests that the Bondi solution fails in this regime and 
the gas flow becomes time-dependent. Indeed, when 
rB = GMh/c^s,oo > '''ta the halo circular velocity is larger 

than the gas sound speed, {GMh/rtaY^"^ = Vdr > Cs,oo- 
Thus, the Bondi solutions transition to "cold" cosmolog- 
ical gas infall, for which self-similar solutions exist. In 
addition, since Hubble velocity at the turnaround radius 
equals the halo circular velocity, rtaH ^ Vdr, we also 
have tH < Tta/ Cs,ao < tcr- The assumption of steady 
flow is therefore no longer correct and the flow becomes 
self-similar and time-dependent. 

Self-similar infall solutions for a mixture of gas and 
dark matter (assuming that the gas is a small fraction 
of the dark matter) are almost identical to the case of 
cold gas infall. The main difference is that well be- 
haved solutions exist for 7 < 4/3 because the gravita- 
tional instability, leading to collapse for 7 < 4/3, oc- 
curs only if the gas is self-gravitating. Two types of self- 
similar solutions exist. If the central black hole absorbs 
all the infalling gas, the solution has free-fall boundary 
conditions at the center (see discussion in § [4]). Most 
of the halo gas is accreted onto the central black hole, 
hence the accretion rate in units of the Eddington rate is 
m = M/MEd ~ tEd/tn- Since tEd ^ 10'' yrs, m » 1 at 
any redshift z }^ 30. The second type of solution is char- 
acterized by shocked accretion of collisional gas. The 
boundary conditions are v = at the center (i.e., no 
accretion onto the central black hole) and pressureless 
gas upstream of the shock. A standing shock occurs at 
a fixed fraction of the turnaround radius rs/rta, whose 
value depends on the polytropic index 7. For 7 = 5/3 
the shock is at — 0.339rtQ, analogous to the result for 
collisionless matter. For smaller values of 7 the shocked 
core is smaller. Self-similar solutions can be found also 
for the cases in which the central black hole accretes an 
arbitrary fraction of the infalling gas. In this case the 
size of the shocked core is further reduced. i Tsuribe et al.l 
(|1995l ) found a self-similar solution for spherical accretion 
around a point mass including radiation drag. They con- 
clude that radiation drag (i) slows down the infalling gas 
velocity far from the center, (ii) flattens the density pro- 
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file near the center, and (iii) reduces the mass- accretion 
rate. The reduction is about a factor of two if the gas 
becomes fully ionized by redshift z ~ 400 — 200. 

Taking into account all the possible accretion 
regimes described above may be quite complex (iLoebl 



19931 iFukue fc Umemural 11994 iMineshige et all Il998l: 
Ciotti et al.l 12004 IWang et all 120061: iPark h Ostriked 
2007() . Fortunately, it is not necessary to estimate in 



detail the accretion rate for the cases in which the flow is 
time-dependent. This is because it is easy to show that 
the accretion rate is always greater than the Eddington 
rate. Feedback effects are likely to regulate the accretion 
rate in this regime. Observations of qu asars at moderate 
redsh ifts and theoretical arguments (jCiotti &: Ostrike"i1 
|2001D suggest that when the gas supply is large (to > 
0.1 — 1), massive black holes shine at nearly the Edding- 
ton rate during a small fraction ~ 1 — 5% of their life, 
while during the rest of the time they are quiescent. Most 
of the energy is emitted during the active phases. Theo- 
retical modeling of feedback processes that produce the 
duty cycle is beyond the scope of the present paper. We 



can treat this case phenomenologically assuming a fidu- 
cial value of the duty cycle inferred from observations. 

In order to demonstrate that when tcr > tn we have 
m ^ 1, it is convenient to rewrite the Bondi accre- 
tion formula in equation ()29|) . in terms of the crossing 
time and Hubble time. Using the relationship fjj = 

{Vldml^b)Gpgas ~ 5Gpgas, and assuming tcr ~ tn and 
A 1, we have 



m ' 



A 



tEdtcr 



» 1, 



at any redshift z > 30. Thus, as anticipated, when 
the accretion flow transitions to a time-dependent self- 
similar regime, the accretion rate is super-Eddington. 
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